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AERONAUTIC SYMBOLS 


1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

. 

Abbrevia- 

tion 

Unit 

Abbreviation 

Length 

Timp 

1 

t 

F 

meter : 

m 

s 

kg 

foot (or mile) 

second for hour) 

ft (or mi) 
sec (or hr) 
Ib 

Force 

weight of 1 kilogram.- 

weight of 1 pound 


P 

V 

horsepower (metric) 


horsepower 

fps 

Speed : 

/kilometers per hour 

(meters per second 

. 

kph 

mps 

miles per hour 

feet per second 


2. GENERAL SYMBOLS 


Weight=»iSf 

Standard acceleration of gravity— 9.80665 m/s* 
or 32.1740 ft/sec* 

W 

Mass= — 

g 

Moment of inertia=mP. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


V Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m~*-s* at 15® C 
and 760 mm; or 0.002378 lb-ft“* sec* 

Specific weight of “standard” air, 1.2255 kg/m* or 
0.07651 Ib/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 

True air speed 

Dynamic pressure, 2 

Lift, absolute coefficient 
Drag, absolute coefficient 


Do 


Profile drag, absolute coefficient 
Induced drag, absolute coefficient Gof=^ 

Parasite drag, absolute coefficient Gdj,=^ 

■' C 

Cross-wind force, absolute coefficient Gc—-^ 


ill. 

it 

Q 

Q 

E 


a 

€ 

«0 

tXi 

OCa 

y 


Angle of settingj)f wings (relative to thrust line) 
Angle of stabilizer setting (relative to thruet 


line) 


Resultant moment , 

Resultant angular velocity 

VI 

Reynolds number, p — where I is a linear dimen- 
sion (e.g., for an airfoil of 1.0 ft chord, 100 
mph, standard pressure at 15° C, the corre- 
sponding^ Reynolds number is 935,400; or for 
an airfoil of 1.0 m chord, 100 mps, the corre- 
sponding Reynolds number is 6,865,000) 
Angle of attack 
Angle of downwash 
Angle of attack, infinite aspect ratio 
Angle of attack, induced 

Angle of attack, absolute (measmed from zero- 
lift position) 

Flight-path angle 
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ANALYSIS OF PLANE-PLASTIC-STRESS PROBLEMS WITH AXIAL 
SYMMETRY IN STRAIN-HARDENING RANGE ^ 

By M. H. Lee Wtr 


SUMMARY 

A simple method is developed jor solving plane-plastic-stress 
problems with axial symmetry in the strain-hardening range 
which is based on the deformation theory oj plasticity employing 
the finite-strain concept. The equations defining the problems 
are first reduced to two simultaneous nonlinear differential 
equations involving two dependent variables: (a) the octahedral 
shear strain, and (6) a parameter indicating the ratio of principal 
stresses. By multiplying the load and dividing the radius by 
an arbitrary constant, it is possible to solve these problems 
without iteration for any value of the modified load. The con- 
stant is determined by the boundary condition. 

This method is applied to a circular membrane under pressure, 
a rotating disk without and with a central hole, and an infinite 
plate with a circular hole. Two materials. Inconel X and 
16-25-6, the octahedral shear stress-strain relations of which do 
not follow the power law, are used. Distributions of octahedral 
shear strain, as well as of principal stresses and strains, are 
obtained. These results are compared with the results of the 
same problems in the elastic range. The variation of load with 
maximum octahedral shear strain of the member is also 
investigated. 

The following conclusions can be drawn: 

1. Inasmuch as the ratios of the principal stresses remain 
essentially constant during loading for the materials considered, 
the deformation theory is applicable to this group of problems. 

2. In plastic deformation, the distributions of the principal 
strains and of the octahedral shear strain are less uniform than 
in the elastic range, although the distributions of the principal 
stresses were more uniform. The stress-concentration factor 
around the hole is reduced with plastic deformation, but a high 
strain-concentration factor occurs. 

3. For the rotating disk and the infinite plate the deformation 
that can be sustained by the member before failure depends 
mainly on the maximum octahedral shear strain of the material. 

4 . The added load that the member could sustain between the 
onset of yielding and failure depended mainly on the octahedral 
shear stress-strain relations of the material. 

INTRODUCTION 

In the design of tiu-bine rotors, it is desirable to know the 
detailed stress and strain distributions in the strain-hardening 


range and the increase in load that can be sustained between 
the onset of yielding and failure. It is also desirable to Imow 
the effects of a notch or a bole in a turbine rotor or other 
machine members that are stressed in the strain-hardening 
range. If a member is thin, it can be analyzed on the basis 
of plane stress. For problems of this type for ideally plastic 
material, Nadai obtained solutions for a thin plate with a 
hole and a flat ring radially stressed (reference 1), and 
Nadai and Donnell obtained a solution for a rotating disk 
(reference 2). For materials having strain-hardening 
characteristics, a solution of plane-stress problems has been 
obtained by Gleyzal for a circular membrane under pressme 
(reference 3). The concept of infinitesimal strain was used 
and the solution was obtained by an iterative procedure 
with a good first approximate solution. The plastic laws 
were always satisfied by using a chart given in reference 3. 
In reference 4, a trial-and-error method is given for a rotating 
disk with very small plastic strain, in which the elastic 
stresses and strains are used as the first approximate values. 
An experimental investigation of high-speed rotating disks 
is given in reference 5 ; distributions of plastic strains (loga- 
rithmic strains) for different types of disk are measured. 
Reference 6 gives an experimental investigation of the burst 
characteristics of rotating disks; stress at the center of the 
disk is calculated by assuming that the material behaves 
elastically at the burst speed; the average tangential stress 
along the radius at burst speed is also calculated. 

A simple method of solving plane-plastic-stress problems 
with axial symmetry in the strain-hardening range for finite 
strains was developed at the NACA Lewis laboratory dming 
1949-50. This method is based on the deformation theory 
of Hencky and Nadai (references 7 to 9), which is derived 
for the condition of constant directions and ratios of the 
principal stresses dming loading. The equations of equilib- 
rium, strain, and plastic law are reduced to two simultane- 
ous nonlinear differential equations involving thi’ee variables, 
one independent and two dependent, that can be integrated 
numerically to any desired accuracy. These variables are 
the proportionate radial distance, the octahedral shear 
strain, and a parameter a that indicates the ratio of principal 
stresses. The magnitude of variation in calculated values 
of the parameter a with change in load directly indicates 
whether the deformation theory is applicable to the problem. 


' Supersedes NAOA TN 2217, “Analysis ol Plane-Stress Problems With Axial Symmetry in Strain-Hardening Range” by M. H. Lee Wu, 1950. 
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The method developed is applied to: (1) a circular 
membrane under pressure, in order to compare results 
obtained by this method with those obtained by Gleyzal 
(reference 3); (2) rotating disks without and with a circular 
central hole, in order to investigate plastic deformation in 
such disks and the effects of the hole; and (3) an infinite 
plate with a circular hole or a flat ring radially stressed, in 
order to investigate the effects of the hole in the strain- 
hardening range. 

In the investigation of (2) and (3), two materials, Inconel 
X and 16-25-6, with different strain-hardening character- 
istics were used in order to determine the effect of the 
octahedral shear stress-strain curve on plastic deformation. 
The octahedral shear stress of these two materials is not a 
power function of the octahedral shear strain, so that more 
general information can be obtained. Distributions of 
stresses and strains of the same problems in the elastic range 
are also calculated for purposes of comparison. 

Acknowledgment is made to Professor D. C. Drucker for 
his discussion of this work and for his suggestion to examine 
whether the logarithmic strain could be applied correctly 
to the present problems and his suggestion to plot the 
stress-strain curves of Inconel X and 16-25-6 on a logarithmic 
scale in order to show that these mateiials do not obey the 
power law. 

SYMBOLS 


The following symbols are used in this report: 

A, B, C, D, E, F coefficients of nonlinear differential equa- 

T 

tions; functions of a, 7, and ^ 


a 

h 

c 

0, H, J, L 
h 

hinit 

k 

V 

r 

s 

u 

w 

z 

a 

7 

e 


e 


initial radius of hole 

initial outside radius of membrane, rotat- 
ing disk, or flat ring 

initial outside radius of plate, very large 
compared with radius a 
trigonometric functions of a 
instantaneous thickness of membrane, 
rotating disk, or plate 
initial thickness of membrane, disk, or plate 
arbitrary loading constants 
constant having a dimension of length 
pressure on membrane 

radial coordinate of undeformed membrane, 
disk, or plate 
arc length 
radial displacement 
axial displacement 
axial coordinate 

parameter indicating ratio of principal 
stresses 

octahedral shear strain 
logarithmic strain (natural strain), loga- 
rithm of instantaneous length divided by 
initial length of element 
angular coordinate 


P 

<7 

T 

CO 

Subscripts; 

b 

c 

0 


1, 2, 3 
r, e, z 


mass per unit volume 

true normal stress, normal force per unit 
instantaneous area 
octahedral shear stress 
angular velocity 

at radius b 
at radius c 

at center for member without hole; at 
radius a for member with concentric 
circular hole 

principal directions in general 
principal directions: radial, tangential, and 
axial directions 


STRESS-STRAIN RELATIONS IN PLASTIC DEFORMATION 

The deformation theory of plasticity for ideally plastic 
materials was developed by Hencky from the theory of 
Saint Venant-Levy-Mises for the cases in which the directions 
and the ratios of principal stresses remain constant during 
loading (reference 7). Nadai extended the theory to include 
materials having strain-hardening characteristics (references 
8 and 9). The conditions for the deformation theory have 
been emphasized by Nadai (reference 9, p. 209), Ilyushin 
(references 10 and 11), Prager (reference 12), and Drucker 
(reference 13). Experiments conducted by Davis (reference 
14), Osgood (reference 15), and others on thin tubes subjected 
to combined loads with the directions and the ratios of the 
principal stresses constant throughout the tube and remaining 
constant during loading show that good results can be 
expected from the deformation theory. 

In more recent experiments on thin tubes by Fraenkel 
(reference 16) and Davis and Parker (reference 17), it has 
been shown that even with considerable variation of the 
ratios of principal stresses during loading the strains obtained 
from the experiments were in good agreement with the strains 
predicted by use of the deformation theory. Further ex- 
perimental investigation is needed to determine the extent 
to which the variation of ratios of principal stresses is 
permissible with the deformation theory. However, when 
the variation is small (approximately 10 percent over the 
strain-hardening range), the deformation theory can be 
expected to give good results. 

In the present problems with axial symmetry, the direc- 
tions of the axes of the principal stresses remain fixed during 
loading and it is probable that the ratios of principal strains 
and of principal stresses also remain approximately constant. 
The deformation theory previously discussed is therefore 
used. The stress-strain relations are as follows: 


ei + «2+es — 0 

(7 X <T 2 O' 2 — <7 3 <7 3 (7 1 

€1 €2 62 63 €3 €1 

T = r(7) 


( 1 ) 

( 2 ) 

(3) 
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where 

r = K 0-1 — 0 - 2 )^ -t- ( q -2 — 0 - 3 )^ + ( 0-3 — (4 a) 

7=|[(ei-e2)^-f(62-.63)^ + (e3-€i)^*'^ (4b) 

From equations (1), (2), and (4a) or (4b), the following 
relations are obtained: 


7['^2-|(<^3+<r,)] 
7[<^3-i(<ri + <72)] 


For plane-stress problems < 73 = 0 . It is convenient to use 
cylindrical coordinates for the problems considered; the 
principal directions 1, 2, and 3 in the preceding equations 
become radial, circumferential, and axial directions, re- 
spectively. The equations thus become 


6r+e(i-|-ej — 0 


and 


T = ^(<r/-<r,<r« + <r«^« 

7 = 2 ..y/.^(er^+CrE9+e9^)^‘ 

e, = i I 


(la) 

(5a) 

(5b) 

(6a) 

(6b) 

(6c) 


When (7r and (xe are expressed in terms of and e$, there is 
obtained 


(Xt — 2 — (2er-f e«) 
7 

£79= 2 — (2€9-t-€r) 

7 


( 7 ) 


to the force divided by the instantaneous area and the 
strains are defined by the following equation: 


5(e^) = 


m 

h 


where Ij is the instantaneous length of a small element having 
the original length of (lj)o and j is any principal direction. 
During plastic deformation, the plastic strains at a particular 
state depend on the path by which that state is reached. 
For the paths along which the ratios of principal stresses 
remain constant during loading, however, the octahedral 
shear stress-strain relation, the value of the octahedral shear 
strain, and the values of the principal strains are defined by 
the initial and final states (references 14, 15, and reference 9, 
p. 209); 5(e^) is then an exact differential and 


It should be noted that the condition under which equation 
(8) was obtained is also one of th^ conditions under which 
the deformation theory is derived; as long as the defonnation 
theory is applicable, equation (8) can also be used. 

EQUATIONS OF EQUILIBRIUM AND STRAINS 
INVOLVING DISPLACEMENTS 
CIRCULAR MEMBRANE UNDER PRESSURE 

Equations of equilibrium and equations of strain are de- 
rived for a circular membrane under pressure. The mem- 
brane considered is so thin that bending stress can be 
neglected (reference 18, p. 576). Figure 1 shows the mem- 
brane clamped at the rim and subjected to a pressure p and 


Og 



Figure 1. — Thin circular membrane (under pressure) and its element in deformed state. 


Because large deformations in the strain-hardening range 
wifi, be considered, the concept that the change of dimension 
of an element is infinitesimal compared with the original 
dimension of the element is not accurate enough. Hence, 
the finite-strain concept, which considers the instantaneous 
dimension of the element, is used. (The equations of in- 
finitesimal strains considered as special cases of finite 
strains are given in appendix A.) The stress is then equal 


a small element defined by A9 and As taken at radius r+u 
in the deformed state. In the undeformed state, the same 
element would be at radius r and defined by A0 and Ar. The 
dotted lines represent an undeformed membrane. The in- 
stantaneous thiclmess of the element and the stresses acting 
on the element are also shown in the figure. Two principal 
stresses are cr, and <xe, and (p is the angle between (x, and the 
original radial direction. 
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Equations of equilibrium. — When all the forces acting on 
the element in the direction of Or are summed up, the follow- 
ing equation of equilibrium is obtained; » 

A0(A-l-A/i,) cos A^-h 


2o-e As(^h+^Ah^sm^ cos <p—pAs(r+u)Ae sin^=0 


When A{r-^u) approaches zero as a lunit, the differential 
equation of equilibrium may be obtained : 


{r+u) 


d{crh) 

d(r+u) 


h,{<n—(rr) 


( 9 ) 


A cap of the membrane bounded by radius r-\-u and the 
forces acting on it are shown in figure 2. Summing up the 
forces in the 2 -direction yields 


taken at radius r+u, defined by Ad and A{r-\-u), are shown in 
figure 3 with aU the external forces acting on the element. 



pTr{r-\-uY=crr 2Th(r+u) 



X 



Figuee 2.— Cap of membrane with radius r + u in deformed state. 


Equations of strains. — ^Inasmuch as the element at radius 
r, defined by Ad and Ar in the undeformed state, is moved by 
the application of pressure p (fig. 1) to radius r+u and 
defined by Ad and As, by use of equation (8) the strains are 


Then 


€r = loge 


dr 


6 9 = log., 


r-\-u 

r 


^=log, 


A 

hinit 


e^r 


d{r-\-u) (. . r dw 
dr ( jd{r-\-v)_ ) 






(11a) 

(lib) 

(11c) 


ROTATING DISK 

Equation of equilibrium. — A disk of radius h and thickness 
h, rotating about its axis with angular speed co, and an element 


Summing up all forces acting on the element in the radial 
direction yields 

(Tf{r-\-u)hAd — (o’r-i-A(rr)[r-l-'it“bA(r-l-u)]A0(/!,-j-AA,)-|- 
2cre[A(r-l-u)] ~ 

pT > , 1 I px[(r-l-Ar)2— r2]A0j 


When A{r-\-u) approaches zero as a limit, the following 
equation of equilibrium is obtained: 


(r+u) 




r-\-u dr 


'd(f+v,)~^'''’ r d(r+u) 

Equations of strains.- — ^The strains are 


( 12 ) 


^ d{r+u) 
dr 

(13a) 

, r-\ru 
e^e — 

r 

(13b) 

^init 

(13o) 


INFINITE PLATE WITH CIRCULAR HOLE OR FLAT RING RADIALLY 
STRESSED 

An infinite plate uniformly stressed in its plane in all 
directions and having a circular hole is shown in figure 4. 
The whole system is equivalent to a very large circular plate 
of radius c with a small concentric circular hole radially sub- 
jected to the same uniform stress tr on the outer boundary. 
The solution obtained for such a plate within any radius 6 
can also be considered as a solution of a flat ring with outer 
radius 6 and inner radius a, that is, uniformly loaded at the 
outer boundary with the radial stress vs obtained in the 
plate solu tion. 

The equations for this case can be obtained in a manner 
similar to the two previous cases, or by simply setting dwjdr 
and w equal to zero in equations for the membrane, or by 
setting 03 equal to zero in the equation for the rotating disk. 
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EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY IN 
TERMS OF PRINCIPAL STRESSES AND STRAINS 

CIRCULAR MEMBRANE UNDER PRESSURE 

A set of ten independent equations (equations (la), ( 3 ), 
(5b), ( 6 a), ( 6 b), (9), (10), (11a), (11b), and (11c)) involving 
the ten unlmowns (t„ ce, «r, e^, y, t, h, u, and w define the 

problem of the circular menlbrane under pressure. If equa- 
tion (lib) is differentiated with respect to r and combined 
with equation ( 11 a), 

des 

r dw (14) 

\ "^L<^(?’+m)J ; 


Substituting equation ( 10 ) in equation (14) to eliminate w 
yields the following equation of compatibility: 


r 






(15) 


Equations (9) and (15) can be simplified by using equations 
( 11 ) to eliminate u and h, which results in 



r~ + arr^ = {(T0—a-,)e^<^r ^ 1 — 


and 


Vp6 ^ 




dr 


L ^^imtO'r _J j 


(16) 

(17) 


The ten equations defining this problem are now reduced to 
seven independent equations, (la), ( 6 a), ( 6 b), (5b), (3), (16), 
and (17), with the seven unloiowns Vr, <re, ee, r, and 7 . 

The solution of the problem is simplified by introducing an 
arbitrary constant k into equations (16) and (17): 




(0-9— a-,)e">- 




h h 


L 2cr, 


r dte 



'PKLe^u->y 

h 

lOinit 

L 2<r, J 


-I2\ 1/2 


— 1 


1/2 Y 


(18) 


where k is any arbitrary unknown constant with the dimen- 
sion of length. By use of the two parameters rjk and 
pklhinit, it is possible to solve the problem in a simple, direct 
way without the use of the iteration. This will be fm'ther 
discussed in the section “Methods of Numerical Integration.” 


ROTATING DISK 


For the rotating disk there are nine independent relations 
(equations (la), (3), (5b), ( 6 a), ( 6 b), (12), (13a), (13b), and 
(13c)) with the nine unknowns <r„ <re, (r, ee, e^, y, r, h, and u. 
If equation (13b) is differentiated with respect to r and 
combined with equation (13a), the following compatibility 
equation is obtained : 


r 


dee 

dr 






(a) Infinite plate with circular hole uniformly stressed in its plane in all directions. 

(b) Flat ring radially stressed. 

(c) Element. 

Figure 4. — Infinite plate with circular hole, flat ring radially stressed, and its element in 

deformed state. 


As in the case of the membrane, u and h can be eliminated 
from the equilibrium equation (12) by using equations (13), 
which yields 

r ^+<rrr ^={<Je-<rr)e^^'-^e)-po,W~^^'> ( 20 ) 


The nine equations defining this problem are now reduced 
to seven independent equations, (la), ( 6 a), ( 6 b), (5b), (3), 
(19), and (20), with the seven unknowns o-r, a-e, e„ ee, e.^, r, 
and y. 

The solution of the problem is made simpler by introducing 
an arbitrary constant k into equations (19) and (20) : 



r dt 




© 


-{<je — <7,)e^*>' — p(o}ky ‘d 


r dee 


'‘da 




(19) 


( 21 ) 
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By use of the parameters rjk and o>k instead of r and co, a 
simple direct solution is possible for any arbitrary value of 
mk with k to be determined by the boundary condition. 

INFINITE PLATE WITH CIRCULAR HOLE OR 
FLAT RING RADIALLY STRESSED 

The equations of equilibrium and compatibility for the 
infinite plate w'ith a circular hole or the flat ring radially 
stressed are: 


(reference 1, p. 189) and for a rotating disk (reference 2). 
From equations (6a), (6b), and (23), the principal strains 
also can be expressed in terms of y and a : 

6r=— ^(sin a— Vs cos a) 

2 V 2 

- (24) 

T • ! — 

69= — ^(sm q;+ V 3 cos a) 

2V2 J 


r 

a 





(<re — <Tr) 





1 




(22) 


The problem is defined by equations (22) together with 
equations (la), (6a), (6b), (5b), and (3) (seven equations 
with seven unknowns). 


EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY IN 
TERMS OF OCTAHEDRAL SHEAR STR.MN AND PARAM- 
ETER INDICATING RATIO OF PRINCIPAL STRESSES 

In the preceding section, displacements are eliminated 
from the equations, which result in seven equations involving 
the seven unknown quantities cr,, as, ee, r, and y. The 
quantity can be expressed in terms of and ej (from 
equation (la)). Two of the fom unlaiowns ar, ae, tr, and e* 
can be eliminated by using equations (6a) and (6b) or (7). The 
quantity t is a known fimction of y that is experimentally 
determined. The problem is then reduced to one involving 
three unlmowns. Obtaining the solution of the resulting 
equations is not, however, a simple matter; the iterative 
process is usually needed. 

It is proposed that this can be avoided by using the fol- 
lowing transformation; 

o-9-t-(Tr=3 V2r sin a 


or 


ai-^ar—^T cos a 
<7r=-Y^T(V3 sin a— cos a) 
ae =-^ sin a-l-cos a) 


(23) 


Then <r, and ae satisfy equation (5a), because the yielding 
surfaces are ellipses according to the deformation theory. 
The octahedral shear stress t, a function of y, in the preced- 
ing equations varies with rjk and also with loading. Such a 
transformation has been used for ideally plastic material 
(T=constant) by Nadai in the section “Yielding in Thin 
Plate With Circular Hole or Flat Kings Radially Stressed” 


The equations of equilibrium and compatibility for the 
three problems considered herein are then obtained in terms 
of y and a in the following form: 


dc 




-C 


doc , dy jp 


(25) 


where the coefficients A, B, C, D, E, and F are functions of 
a, y, and rjk. For the circular membrane under pressure, 
from equation (18), 

H=(V3 cos a-psin a)— (V3 sin a— cos a) 

V2 


K=(-\/3 sin a— cos a) 


y dr y sin a\ 1 

V2 / 


(7=2 (cos a) e 


( ^llrcosa') 


„Vf(V3 


sin a-}-C 0 S a)y 


L 6r^(V3 sin a— cos a)‘ 
Z1 = (V3 sin a — cos a)y 
E=—{-\j2> cos a-bsin d) 

F=2V2 


/ r V/ pk 

\k) {kini, 


'pk'^j 


1/2 


e 


(-Videos a) 


.VI (V3 


sin CK+COS a) 7 


L 6r^(V3sina; — cos a)^ 

For the rotating disk, from equation (21), 




j 


l/ 2 \ 


-• ' J 


(26) 


y cos a 


H=(V3 cos a-bsin a)—{-yjZ sin a— cos a) ^ 
^y dr y sin a 


r dy ^ 
p{oikf 


sin a— cos d) 

C=2(cosa)e^ V§'>'‘=°®“)_^2 

i)=(V3 sin a— cos a)y 
H=— (V3 cos a-bsin a) 

F=2V2[i-/ Vr“'“)] 




1 /?-V_ 


y (27) 
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(28) 


For the infinite plate with a circular hole, from equation (22), 

-4= (Vs cos a-fsin a) — (VS sin a— cos a) “ 

-y 2 

, I— , \ /y dr j sin a\ 1 

S= (VS sin a-cos ^)y 

C=2 (cos a) 

Z>=(VS sin a— cos a) y 
E=-(-)jd cos a-l-sin a) 

F= 2 V 2 [l-e("^"“^“)] 

With these transformations, the solution of the problems 
is reduced to simply a numerical integration of the two simul- 
taneous differential equations (equations (25)) involving the 
two unlmowns y and a. Furthermore, the parameter 7 , 
being the octahedral shear strain, directly indicates the 
stage of plastic deformation at any point under any load. 
(In plastic problems, according to the deformation theory, 
the individual stress and strain distributions cannot give as 
clear a picture of the stage of plastic deformation as can the 
octahedral shear strain.) Also, the parameter a indicates 
the ratio of the principal stresses or strains. At any point, 
if a remains constant during loading, the ratio of principal 
stresses at that point remains fixed. The value of a obtained 
at each point in the calculation during loading directly indi- 
cates whether or not the deformation theory is applicable to 
the problem. 

The value of a is known at the boundaries or the center. 
This value can be determined from equations (23) and (24). 
For the circular membrane under pressure, 
when r /6 = 0 , 

<7 


when rlb = l, 


a=^=1.5708 


a=^ 7 t=2.0944 


For the rotating disk without a hole, 
when rlb = 0, 

(77-= <70 


when r /6 = l, 


a=V=1.5708 


(7, = 0 


o:=;^=0.5236 

6 


For the rotating disk with a hole, 
when r/a=l and rlb = l, 


(Tr=0 

01 = 1 = 0.5236 

6 

For the infinite plate with a circular hole, 
when r!a=\, 

CTr = 0 

a=|=0.5236 

D 

when rja approaches cja or a value that is large compared 
with 1 , 

(T j — (X0 

a=|=1.5708 

METHODS OF NUMERICAL INTEGRATION 

Two methods are developed to solve the differential equa- 
tions (25). In the first method, the differential equations 
are numerically integrated along rjk, which is considered the 
independent variable. (In the second method, ot is con- 
sidered the independent variable.) Because many terms in 
the equations are trigonometric functions of a, the use of a 
as the independent variable considerably reduces the work 
of computation. 

Numerical irtegration with r/fc as independent variable. ^ — 
Equation (25) can be written in the following forms: 


r do( 

CE-FB 


AE-DB 

r dy 

FA-CD 


EA-BD 


At any point, if a and 7 are known. 


and 


dy 


dirjk) dirjk) 


(29) 


can be 


calculated by equations (29). At the boundaries or the 
center, a. is luiown, but 7 is determined by the load. Only 
one value (unlcnown) of 7 corresponding to a particular load 
exists on each boundary; therefore it is difficult to start the 
numerical integrations on the boimdary with the correct 
value of 7 corresponding to a given load. Also, in plastic 
problems covering the strain-hardening range, the method 
of superposition is invalid. Usually, a method of iteration 
is used to solve the problem (for example, references 3 and 4). 
In the method presented herein, an arbitrary but unlcnown 
constant k has been introduced in equations (18), ( 21 ), and 
( 22 ). For the cases considered, the terms in the equations 

that involve load are always multiplied by r, so that ^ 
can be written as ^ ^ in equations (18) and (26) 

and {(ary as {<aky in equations (21) and (27). 


969882-51 2 
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The numerical integration then can be started at the 
inner boundary (or at the center if there is no circular hole 
at the center) by using the known values of a desired 


value of To, 


and an arbitrary value 



for the 


membrane or of (uky for the rotating disk. The munerical 
integrations can then be carried out, obtaining values of a 
and y at different values of rjk, until a progressively reaches 
the value that satisfies the outer boundary condition. 
Because the value of r is Imown at the boundaries, the value 
of k can be determined for the selected value of yo- The 
number of points and the formulas used in the calculation 
depend on the accuracy required (references 19 and 20). 
If the formula for evaluating definite integrals is applied 
after using the forward integration formula (references 19 
and 20), great accurac}^ can easily be obtained. 

The procedure used herein to obtain solutions is the same 
for each problem. Calculations are started from the inner 
boundary (or from the center if there is no circular hole at the 
center) with the known value of a„, the desired value of 
7„, and the arbitrary loading term. The parameter <*„ is 
equal to tt/2 at r/6 = 0 for the membrane and for the solid 
rotating disk and is equal to x/6 at r/a=l for the infinite 
plate with a circular hole and for the rotating disk with a 


hole. 


The arbitrary loading terms are 



and {wky for 


the membrane and the rotating disk, respectively. Then 

corresponding to and y„ at the 

inner boundary or the center, are obtained from equations 
(29). The following formulas for forward integration are 
used to determine the first approximate values of a and y 
at the next point (ui* and 7i*) ; 


t doi ~\ , r ^7 ” 

d{rlk)i Ld(r/fc)J„ 


again in order to calculate the values of 


r dc 


r dy 

Ki)J, 


K0J. 


and 


. By use of the following formulas for evaluating 


definite integrals, the volumes of ai and 71 are calculated: 


1 1 

“ 1 
s-l -j 

■a 

1 1 

da 

l<i)J 

+ 

0 

da 

ki)J 

] 

[(§).-(a 

dy ~ 

+ 

0 

~ dy ~ 

K0J 

]| 


y (30c) 


This procedure is applied to the next point, and so forth, 
until the value of a reaches the required value of a* at the 
outside boundary (0:6= 2/3 v at r/6 = l for the membrane, 
0:5= x/6 at r/6 = l for the rotating disk, and qi(,= x /2 at 
rja=cja for the thin plate with a circular hole). Inasmuch as 

/r\ 

\kj'^=‘‘i, k 

the loading terms are determined as follows: 

For the membrane, 



(31a) 


For the rotating disk. 


(o>by= 



(31b) 


«i* = ao + |^ 

.(0, (0.- 

da 

<0 

0 

Ti* = 7o+|^ 


1 1 

73 ^ 

J 1 1 

> 

", 


(30a) 


By substitution of a^* and 71* into equation (29), approxi- 
da “1 j r dy 


mate values of 


and 


J, K0J. 


are obtained and 


the second approximate values of a, and 71 (a/ 
can be computed from the following formulas: 


and 


7i 



The values of «i** and 71** are substituted into equation (29) 


For the infinite plate with a circular hole, 

Te 

tc=crhc (31c) 

or for the flat ring radially stressed at the outside diameter 

b, 


tjj — (o'r)&^& 


(^)r< 


= i<^T)bhinue 


(sin cos a^) 


(31d) 


where tc- and tb are the tensions per unit original circum- 
ferential length at r=c and r = b, respectively. 

Numerical integration with a as independent variable. — 
Equations (29) can be written in the following forms: 


dy FA-CD ' 
da~ CE-FB 

AE-DBr 
da CE—BF k^ 


( 32 ) 


By use of equations (26) to (28) and expansion of 
into a series, the following equations are obtained: 
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For th,e circular membrane, from equations (26), 


CE-BF== 26 J-\^ -^ZHJg (a, ^)/i («, 7)+2Vs] 
AF-CD ^2 42 L- 2HJj -242Lf^ (a, t) j («, 7, ^ 
AE-BD=-n-J^g(^a,y,^^) 


./ r pk\ ^ 


(33) 


For the rotating disk, from equations (27), 

( 7 F;- 5 F= -2Hi- 2 V3HJ3 («, -y, J (a, 7) + i^(0 /s (a, 7) 

AF-<7£>= J 8H^- 2 V3i?i [1 (a, 7 )] + («, 7 ) |7 

AE~BD=-L>-J>s(o:,yi<^) 

For the infinite plate with a circular hole, from equations (28), 

CE-BF=-2HL-24zHJg (a, t, J ^)/i («, 7) 

AF-CD={8H^-243HL[1-Mc‘,y)]]y 


(34) 


(35) 


where 


6? = sin a , 
i7=cos a 


J - 43 sin a— cos a 
L= 4 ^ cos a+sin a 



and 

/i(“,7) = 






- (cos a)y 



V 


^ (cos a)7+;j (cos^ a)y^ 


= 1 — -y/| (cos a)yfi{a, y) 

— ~ sin a 

Uo‘,y)=e^^ 

K’i/I™")' 

/ Y d,T\ y dr /3 
./ r pk\ 


y^+ 


7 '+ . . . 


sin a— cos a 

-|t(-/ 3 sin a+cos a) 


_7 d 
T dy 


Qt\4^ si’l a— cos aY 

fi rX M^,y) T^ 

1 ^ ej^KrkJ Ua,y)J 


i-Vi 


7 

J 


J 


The symbols G, H, J, and L are trigonometric functions of 
a only; Ki and are arbitrary loading constants. The 
symbols /i, /j, /s, and g are functions of a. and 7 ; j is a 
function of a, 7 , and rlk. 


r, y -■J^(eOSa)y 

f2{a,y)=e V2 

= 1— -y/l (cos a)7+ix| (cos^ a)y^— 
i “)7®+ ■ • • 


For the solution of an infinite plate with a circular hole a 
is the independent variable. The procedure of numerical 
integration is similar to that used in the method in which 
r/k is the independent variable. The first four terms of the 
series of are used; the accuracy of the result is the same 

as that of the previous method, but computation is reduced 
by one half. 
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Both methods presented herein are used to obtain the solu- 
tions for the given values of y„. The purpose of the present 
paper is to obtain solutions for the entire strain-hardening 
range and the methods developed are very convenient for this 
purpose. If, however, a solution for only a particular value 
of loading is required, it can be obtained by interpolating 
between values obtained from two or three solutions corre- 
sponding to loading near the specified value. 

NUMERICAL EXAMPLES 

Membrane. — In order to compare the results for the 
circular membrane obtained by the method developed herein 
with those obtained by Gleyzal (reference 3), one numerical 
solution for infinitesimal strain is calculated by using the t{j) 
curve of the tensile test in figure 1 of reference 3. Inasmuch 
as reference 3 states that: “For simplicity, strain will be 
taken to mean conventional strain (ds—dso)ldSo where ds and 
dso are final and initial arc length, respectively.”, equations 
(26a) and (36) given in appendix A for infinitesimal strain 
are used. The calculation is started at r/A:=0.005. Values 
of aa=1.6708, 7a=0.0299, and p/:/Ainit= 55,920 are used. 

Rotating disk. — Numerical solutions for the rotating disk 



(a) Linear-scale plot. 



for finite strain (equations (25) and (27)) are calctdated. 
The t(y) curves of two materials. Inconel X and 16-25-6, 
are plotted in figure 5 (a). These data were supplied by 
W. F. Brown, Jr., H. Schwartzbart, and M. H. Jones. The 
same 7 ( 7 ) curves are plotted on logarithmic coordinates 
in figure 5 (b). These materials. Inconel X and 16-25-6, 
for which t is not a power function of 7 , were chosen so that 
more general information can be obtained. The given 
octahedral shear stress-strain curves (fig. 5) of these two 
materials have not been corrected for the triaxiality and 
nonuniform stress distribution introduced by necking and 
consequently do not represent the exact stress-strain relation 
after necking of these two materials. The solutions obtained 
from the t( 7 ) curves of the tensile test after necking can, 
however, represent the solutions corresponding to materials 
having the exact rfy) curves shown in figure 5 and for simplic- 
ity such materials are herein referred to as “Inconel X” 
and “16-25-6”. 

The calculation for the solid rotating disk is started at 
r/jfc=0.005, as in the case of the membrane. 

Three solutions are also obtained for the rotating disk with 
a central hole, using Inconel X. Calculations are started 
at rja=\. 

All numerical examples for the rotating disk are given 
in the following table: 


Solid rotating disk 

Material 

7o 

K2 = -^^p{tak)^ 

Inconel X 

0. 04 

IXIO* 1 


.1152 

1X10« 1 


.30 

1X105 1 

16-25-6 

0. 04 

IX 105 


.1152 

1X105 


.30 

2.5X105 


Rotating disk with central hole 

Material 

78 


Inconel X 

0. 30 

1X10< 


.30 

2X103 


.30 

4X1Q3 


Infinite plate with circular hole.^ — The calculations for 
the infinite plate with a circular hole are carried out for the 
case in wliich cr,=0 at r/a=l. The value of a„ at rja—l 
is then equal to 0.5236. (When is different from 0 at 
r/a=l, the corresponding value of should be used.) The 
same materials as in the previous problem are considered. 
The numerical examples are: 


Material 

7o 

Inconel X 

0. 04 


.1152 


.1871 


.30 

16-25-6 

0. 04 


.1871 


.30 


(b) Logarithmic-scalo plot. 

Figure 5.— Octahedral shear stress-strain curves. 





PLANE-PLASTIC-STEESS PROBLEMS WITH AXIAL SYMMETRY IN STRAIN-HARDENING RANGE 


11 



Figure 6.— Vai'lution ot principal stresses with proportionate radius distance. 


RESULTS AND DISCUSSIONS 

The radial and circumferential stresses dr and ae, respec- 
tively, obtained for the circular membrane are plotted against 
rjh in figure 6. Two curves, taken from reference 3, cor- 
responding to calculations for about the same pressui’e used 
in the present calculation, are included in the figure for com- 
parison. In the present calculation, the t( 7 ) emwe given in 
figure 1 of reference 3 and the same infinitesimal-strain defini- 
tion based on the original dimension are used. The initial 
thickness is also used for consistency in the calculation 
rather than the instantaneous thicloiess h, which is used in 
reference 3. 

The variations of a with the radius for the rotating disk 
and with the radius for the infinite plate with a circular hole 
are plotted in figures 7 (a) and 7 (b) , respectively, for different 
loads and materials. The variations of a with y„ (or loading) 
at various radii for the rotating disk and the infinite plate 
with a circular hole are plotted in figures 8 (a) and 8 (b), 
respectively. Similar curves for the ratio of the principal 
stresses o-r/o-e are shown in figures 9 (a), 9 (b), and 10. Com- 
parison of figure 7 and figures 9 (a) and 9 (b) shows that the 
variations of a with radius are very similar to the variations 
of Crjce with radius, although the relation between a and 
cTr/cre is not linear. 

Numerical examples for a membrane with a large strain are 
not calculated herein, because the result of reference 3 is 
sufficient to give an approximate variation of the ratios of 
principal stresses along the radius during loading, although 
the infinitesimal-strain concept is used. The variations of 
the ratio of principal stresses with radius for different loads, 
based on the values of o-, and ae given in figures 8 and 9 of 
reference 3, are calculated and plotted on figure 9 (c). 

959882—51 3 



(a) Rotating disk. 



(b) Infinite plate with hole. 

Figure 7. — Variations of parameter a \vitfa proportionate radius distance for Inconel X and 

16-25-6. 


•The values of Vr are plotted against <re at various radii 
under different loads for the rotating disk and the infinite 
plate with a circular hole in figures 11 (a) and 11 (b). The 
heavy solid and dashed curves represent the values of <Xr 
and (T* at different radii for any given load and are called 
loading curves. The loading curve moves away from the 
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Maximum ocfahedrai shear sfrain. To 


(a) Rotating disk. 

(b) Infinite plate with hole. 

FiorRE 8.— Variation of parameter a with maximum octahedral shear .strain at different radii. 

origin with increasing load. The light solid and dotted lines 
connectmg the different loading curves at a given radius and 
extending to the origin represent the values of <rr and ae at 
different loads for any given radius and are called loading 
paths. Also shown in the figures are the yielding surfaces, 
which are ellipses under the deformation theory. 

A clear picture of the variation of the ratios of principal 
stresses in this group of problems with different loads for 
Inconel X, 16-25-6, and the material used in reference 3 is 
given in figures 7 to 1 1 . It is evident that the ratios of 
principal stresses remain essentially constant during loading. 



0 J .2 .3 .4 .5 .6 .7 .8 .3 LO 


Proportionate radial distance, rjb 
(ii) Rotating disk. 



Proportionate radial distance, r ja 
(b) Infinite phite with hole. 

I'iGtiRE 9.— Variations of ratio of principal stresses with proportionate radial distance. 


The deformation theory is therefore applicable to this group 
of problems, at least for the materials considered. 

The variations of 7 and 7 / 7 ^ with radius are plotted in 
figures 12 and 13, respectively, for the rotating disk and the 
infinite plate with a circular hole. It is interesting to note 
that the curves in figure 13 for different loads for the same 
material are quite close together. For different materials, 
the curves of figures 7 and 9 are also close, but the curves of 
figure 13 are not as close together. 




Ratio of principal stresses, dpjdg I? Ratio of principal stresses^, dr j(Sg 


7 
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Proporfionofe rodiol distance, r /b 


(c) Circular membrane under pressure. 


JURE 9.— Concluded. 


Variations of ratio of principal stresses with proportionate 


radial distance. 



0 .05 JO J5 .eo .25 .30 


Maximum octahedral shear strain^ 

(a) Rotating disk. (b) Infinite plate with hole. 

Figure 10. — Variation of ratio of principal stresses with marimum octahedral shear strain at 

different radii. 



(a) Rotating disk. 

(b) Infinite plate with circular hole. 

Figure 11. — Loading curves, loading paths, and yielding surfaces. 


The distributions of principal stresses and principal strains 
along the radius for the rotating disk and for the infinite plate 
with a circular hole are plotted in figures 14 and 15, respec- 
tively. For comparison, the variations of ere/ (ve) e, ee/ (ee) t, and 
t/Y 6 with radius for both the elastic and the plastic range are 
plotted in figures 16 and 17. (The equations for the elastic 
case are given in appendix If only the stress distribu- 

tions for the elastic and plastic cases are compared (figs. 

16 (a) and 17 (a)), it is seen that the stresses are more uniform 
in the plastic state; but if the distributions of the principal 
strains and the octahedral shear strain for thb elastic and the 
plastic cases are compared (figs. 16 (b), 16 (c), 17 (b), and 

17 (c)), it is evident that a less-rmiform strain distribution is 
obtained in the plastic state. It is of special interest in the 
case of the finite plate with a hole to note that with plastic 
deformation and stress- (tangential stress) concentration 
factor aroxmd the hole is reduced; instead there is a high 
concentration in principal strain and in octahedral shear 
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(a) Kotating disk. (b) Infinite plate with hole. 

Figure 14.— Variation of principal stresses with proportionate radial distance. 


strain. A similar conclusion regarding the concentration 
factor around a circular hole in a tension panel is given io 
references 21 and 22. 

The quantities (rel(ffs)oi er/(er)o, and ee/(6fl) „ along the 

radius for the rotating disk and (xbI{<ib)o and (bIMo for the 
infinite plate with a circular hole are plotted in figures 18 
and 19, respectively. The curves representing Crl{<Tr)„, erl{fr)o, 
and es/(e«)o for Inconel X and 16-25-6 and different values 
of 7 o are close together; but the curves of cbI{(^b)o are quite 
far apart for the two materials, as well as for different values 
of y„. 

The relation between the rotating-speed function p{t^by 


and 7 „ for the rotating disk and the relation between the 
tension per unit original circumferential length tbihinu and 
7 „ for the infinite plate with a hole are plotted in figures 20 (a) 
and 20 (b), respectively. It is shown in these figures that 
pioihy and increase considerably for Inconel X and 

increase only slightly for 16—25-6 as the value of jo increases 
from 0.04 to 0.30. 

Figures 7, 13, and 16 to 19 show that for the plate with a 
hole, the variations of a, y/jo, erl(er)o, and tBl{eB)o with radius 
are essentially independent of the value of y„ for the plate 
and the T{y) curve of the material, at least for the materials 
considered. These results show that the deformation that 



Principal strains, eg Principal strains 
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proportional radial distance, rjb 


(a) Rotating disk. 



(b) Infinite plate with hole. 

Figure 15.— Variation of principal strains with proportionate radial distance. 



(a) Variation of ratio of tangential stress to minimum tangential stress at outer radius with 

proportionate radial distance. 

(b) Variation of ratio of tangential strain to minimum tangential strain at outer radius with 

proportionate radial distance. 

(c) Variation of ratio of octahedral shear strain to minimum octahedral shear strain at outer 

radius with proportionate radial distance. 

Figure 16. — Comparison of results obtained in elastic and plastic range for rotating disk. 
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(a) Variation of ratio of principal stresses to tangential stress at outer radius with 

proportionate radial distance. 

(b) Variation of ratio of principal strains to tangential strain at outer radius with proportionate 

radial distance. 

(c) Variation of ratio of octahedral shear strain to octahedral shear strain at outer radius with 

proportionate radial distance. 

Figure 17. — Comparison of results obtained in elastic and plastic range for infinite plate with 

circular hole. 


(a) Variation of ratio of principal stresses to stress at center of disk with proportionate radial 

distapee. 

(b) Variation of ratios of principal strains to strain at center of disk with proportionate radial 

distance. 

Figure 18. — Comparison of stress and strain distributions for different materials and for 
different maximum octahedral shear strains of rotating disk. 
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(a) Variation of ratio of tangential stress to tangential stress at inner radius with proportionate 

radial distance. 

(b) Variation of ratio of tangential strain to tangential strain at inner radius with proportionate 

radial distance. 

Figure 19.~Comparison of stress and strain distributions for different materials and for 
different maximum octahedral shear strains of infinite plate with circular hole. 

can be accepted by the plate before failure depends mainly 
on the maximum octahedral shear strain (or ductility) of the 
material, which would not be true if the strain distributions 
were a function of the 7(7) curve. For the rotating disk, 



(a) Solid rotating disk. 

(b) Flat ring radially stressed. 

Figure 20. — Relation between rotating-speed function and maximum octahedral shear strain 
for disk and relation between tension per unit original circumferential length for plate. 

however, a slight effect of 7^ and the 7(7) curve is apparent 
on the strains; this effect seems to be caused by the body- 
force term of the disk. 

The stress distribution that will determine the load which 
a member can sustain is now considered. Figures 16 to 19 
show that the variation of aelMo with radius depends on 
the 7(7) curve of the material and on the value of y„ for the 
member. Figure 20 indicates that the load also depends on 
the 7(7) curve. It therefore follows that the added load 
that the member can sustain between the onset of yielding 
and failure depends on the 7(7) curve of the material. The 
octahedral shear (or effective) stress and strain curve of the 
material should be used as a criterion in selecting a material 
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Ratio of Outer 
■ to in ne r radius 


r -Solid rotating disk 
\ I I I 


SSOxlO^- 


Ratio af outer 
to inner radius 
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■ 6 .8 

Proportionate radio I d istonce, r/b 

(a) Variation of parameter a with proportionate radial distance. 


Ratio of outer 
to inner rodia 


.6 .6 

Proportionate radial distance, rfb 

(b) Variation of octahedral shear strain with proportionate radial distance. 


rnmum 


.6 .a 

Proportionate radial distance, r’fb 

(c) Variation of principal stresses with proportionate radial distance. 


Figure 21. — Rotating disk with hole. Inconel X; -yo. 0.3000, . 









20 


EEPOET 1021 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(d) Variation of principal strains with proportionate radial distance. 

Figure 21. — Continued. Rotating disk with hole. Inconel X: yo , 0.3000. 

for a particular member under a particular loading condition, 
because consideration of the maximum octahedral shear 
strain only (or ductility only) of the material is insufficient. 

The variations of a, y, <Tr, <re, «r, «e, and 7/75 with radius for 
three rotating disks with a hole are shown in figure 21 . The 
values of the ratios of outer to inner radius h\a of these three 
disks equal 6.32, 12.45, and 28.12. These disks were made 
of Inconel X and had a maximum octahedral shear strain 
yo of 0.3 at the inner radius of the disk. The tangential 
stress < 79 , the tangential strain ee, and the octahedral shear 
strain 7 are much less uniform for the disk with a hole than 
for a solid rotating disk. The ratio of maximum to minimum 
octahedral shear strain yohb is equal to 7.41 for a disk with 
bla— 5.32, 11.75 for a disk with 6/a=12.45, and 14.1 for a 
disk with b/a=28.l2; for a solid disk of the same material, 
the ratio 70 / 7 is about 5.3. 

The load, rotating-speed function p{o>hY, for disks of 
Inconel X reachmg a maximum octahedral shear strain y„ of 
0.3 at the inner radius of the disk and having different ratios 
of inner to outer radius ajh is represented by the solid curve 



(e) Variation of ratio of octahedral shear strain to minimum octahedral shear strain at outer 
radius with proportionate radial distance. 

Fioure 21. — Concluded. Rotating disk with hole. Inconel X; 7 - 0 , 0.3000. 



Figure 22,— Variation of load (function of speed) with ratio of inner to outer radius of rotating 
disk with hole. Inconel X; to=0.3000. 
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in figure 22 . The dashed curve in figure 22 is obtained by 
extending this solid curve toward a/ 6 =l, where the value of 
p(caby can be determined by considering a rotating ring with 
a/ 6 ->l. The figure indicates approximately how the load 
p(ci}by varies with disks having different ratios of inner to 
outer radius and reaching the same maximum octahedral 
shear strain, at . the inner radius, of the disk. The value of 
p(o}by for a sohd rotating disk made of Inconel X with 
7o=0.3 at the center of the disk is indicated in the same figure. 

CONCLUSIONS 

The results obtained for a membrane, a rotating disk 
without and with a hole, and an infinite plate with a hole 
strained in the strain-hardening range in which the elastic 
strains are negUgible compared with the plastic strains for 
Inconel X and 16-25-6 in the absence of time and tempera- 
ture effects and unloading show that : 

( 1 ) The method developed not only accurately solves 
the plane-plastic-stress problems with axial symmetry in a 
simple manner but also shows clearly the octahedral shear 
strain distribution and the ratio of principal stresses during 
loading. 

( 2 ) The ratio of the principal stresses in the cases investi- 
gated remained essentially constant during loading and, 
consequently, the deformation theory is applicable to this 
group of problems for the materials considered. 


(3) The distributions of principal strains and octahedral 
shear strains in the plastic state are less uniform than those 
in the elastic state, although the distributions of tangential 
stresses appear more uniform in the plastic state. The 
stress concentration factor around a hole is reduced in the 
plastic state, but instead there is a high concentration of 
principal strain and octahedral shear strain. 

(4) The ratios of the strains along the radius to their 
maximmn value are essentially independent of the value of 
the maximum octahedral shear strain of the member and 
the octahedral shear stress-strain curve of the material. 
Hence, the deformation that can be sustained by the member 
before failure depends mainly on the maximum octahedral 
shear strain (or ductility) of the material. 

(5) The stress distributions depend on the octahedral 
shear stress-strain curve of the material. Hence, the added 
load that the member can sustain between the onset of 
yielding and failure depends mainly upon the octahedral 
shear (or effective) stress-strain curve in the strain -hardening 
range of the material. 

Lewis Flight Peopulsion Laboratory t 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 28, 1950 


APPENDIX A 


EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY FOR 
INFINITESIMAL STRAIN IN TERMS OF a AND y 

The final forms of the equilibrium and compatibility 
equations for small strains are given in this section. The 
concept of infinitesimal strain is defined as follows: The 
changes of dimensions are small compared with the original 
dimensions but are large enough so that the elastic strain 
can be neglected. The equations presented can be obtained 
either by direct derivation as was done previously or by 
reducing the equations for finite strains through expanding 
the terms in series and neglecting the small terms. 

For infinitesimal strain, the coefficients (functions of a and 7 ) 
A, B, C, D, E, and F of equations (25) are each denoted by 
a superscript prime but the coefficients (functions of a and 7 ) 
are simpler than those for large strain. 


^ k 


-B 


dy 


© '"© 


-C' 




I Tpt ^ 

k 


dy 


(26a) 


\k. 

For the circular membrane under pressure, 
A' — ©3 cos a-f-sin a 

B' =(V^ sin a— cos a)i ^ 


C =2 cos a 


D'=(-\f3 sin a— cos a )7 
E' = — (Ys cos a+sin a) 

F' =2 V 37 cos ct+^ 


pk r_ 
hinit k 


6 L'tCVs sin^a — cos^a)_ 


(36) 


For the rotating disk, 

A'=-y[Z cos a-l-sin a 

B'=(J3 sin a — cos a) — 

T dy 

O'-2cos„-.y|pW(0’i _ 

sin a— cos ©7 
— (Vs cos a-fsin a) 

F '=2 Vs (cos ©7 

For the infinite plate with a circular hole, 


(37) 


A'= Vs cos a-j-sin a. 

H'=(VS sin a— cos © — ^ 
T dy 

C' = 2 cos a 
H'=(Vs sin a — cos © 
E'=-{-y[Z COS a-}-sin © 

F' = 2VS(cos©7 


(38) 


For small strains, the coefficients A',B',0', E' ,E' , and F' 
are used in equation (29) instead of A, B, O, D, E, and F, 
respectively. 
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APPENDIX B 


EQUATIONS FOR ROTATING DISK AND INFINITE PLATE 
WITH CIRCULAR HOLE IN ELASTIC RANGE 

KOTATING DISK 

For a solid rotating disk with the radial stress at the 
periphery (r=b) equal to zero, the principal stresses can be 
expressed by the following equations (reference 23, p. 68) : 




(Tg — - 


2,; l+3y 
8 ^ 8 


9 9 


(39) 


where v is Poisson’s ratio. At r—b, 


1 


( 0 - 9 ) 6 =^ P03^b\l — v) 


Dividing equation (39) by (0-9) „ yields 

Cfr 


(<79)6 2(1- 

0-9 3 + 


(0-9)9 2(1 


tlL fi _i±^" (L\"~ 

— (/) L 3 + v \b) _ 




(39a) 


The stress-strain relations of plane-stress problems in the 
elastic range are : 


1 . 

«r = ^ (O-,— 1/0-9 


«9=^ (0-9— J/0-,) 


(40) 


where E is the modulus of elasticity in tension and com- 
pression. 

Substituting equations (39) into equations (40) yields: 




(40a) 


or 


{^i)b 2 

69 3 -|- 

(69)9 2 




(40b) 


The equations for the octahedral shear stress and strain 
given by equations (4a), (4b), and (5a) can be applied to both 
the elastic and the plastic ranges, but equation (5b) can be 
applied only in the plastic range. The octahedral shear 
strain in the elastic range can be calculated by equation (4b) 
or by using the following equation: 


2(1 -fr) _2(1 + j/)V2 


(tr,.^ — o-rO-g-f-o-/)'^^ (41) 


' E E Z 

Substitute equations (39) in equation (41) to obtain: 

_ V2 (l + r ) - 2^2', 

^ 12£'(3 + r) ^ ^ 

[(3 + vf- 4(1 + r)(3 + v) (0'+(7 + 2 r+ 7 r^) (4 1 a) 

or 


76 2(1 — 1 /) 

[(3 + r)^-4( 1 + r)(3 -1- v) (^) V (7 + 2 r+ 7 v^) (0]'"' (41b) 

The value of Poisson’s ratio v for the two materials are: 
r=0.29 for Inconel X (reference 24) 
v= 0.286 for 16-25-6 (reference 25) 

INFINITE PLATE WITH CIRCULAR HOLE 

For a uniformly loaded infinite plate with a circular hole, 
the principal stresses are (reference 23, p. 56); 


<^r==p+2t7 


-A 


0-9 = 


■20 


(42) 


where A and 0 are arbitrary constants. For the plate con- 
sidered herein, tlie boundary conditions are: 

o-, = 0 at r = a 

<Tr=(<Tr)d at 6 

These boundary conditions are used to determine the 
arbitrary constants A and C, which yield 


0-,: 


(V.)9 Q 


1 - 


(I) a 

("£\ 


-fl 


(42a) 


or 


’-(I) (£) J 

lA- 

0-, __ 1 \a/ 

ey 


0-9 


. ey 


+1 


(0-9)9 


1+ 


/ a Y / r Y 

\bj \aj 


(42b) 
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Substituting equations (42a) into equations (40) yields 


or 


1 

©■ 




1 (<7 

€0=t^ 




«r 


© 


(•.). [a~.)+(i+.)(|)-]©- 

d-.)(0+(i+.) 

Substituting equations (42a) into equation (41) yields 


(43) 


(43a) 


or 


2V2(l + r) 

m 





(44) 
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